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I.  Introduction 


The  U.S.  Army  Ballistic  Research  Laboratory  (BRL)  has  been  conducting  research  to 
support  the  design  of  a  facility  intended  to  subject  full  scale  military  equipment  to  the 
blast  and  thermal  environment  produced  by  tactical  nuclear  weapons.  This  facility,  termed 
the  Large  Blast/Thermal  Simulator  (LB/TS)  is  designed  to  produce  a  wide  range  of 
nuclear  blast  and  thermal  environments.  Shock  overpressures  in  the  simulator  can  range 
from  2  —  35  psi  and  simulated  weapon  yield  can  range  from  1  —  600  AfF.  The  design  of  the 
LB/TS  is  being  managed  by  the  Army  Corps  of  Engineers  for  the  Defense  Nuclear  Agency. 
Construction  is  scheduled  to  begin  in  1991. 

The  LB/TS  can  be  thought  of  as  a  large  shock  tube  with  a  complex  geometry.  The  blast 
portion  of  a  nuclear  blast  simulation  is  accomplished  by  releasing  high  pressure  gas  from  a 
set  of  relatively  small  steel  driver  tubes  into  a  large  concrete  expansion  section.  By  adjusting 
the  initial  driver  parameters  of  pressure,  temperature  and  volume,  a  desired  combination  of 
shock  overpressure  and  simulated  weapon  yield  can  be  obtained.  Baffles  are  used  to  adjust 
the  area  ratio  between  the  driver  tubes  and  the  expansion  tunnel,  which  also  affects  the  blast 
simulation.  This  technique  is  employed  to  extend  the  facility’s  operating  range. 

When  completed,  the  driver  system  of  the  LB/TS  will  consist  of  nine  steel  cylinders. 
Each  of  these  cylinders  will  have  an  inside  diameter  of  1.83  m  and  will  have  a  converging 
nozzle  and  diaphragm  system  at  its  downstream  end.  The  flow  in  the  simulator  is  initiated 
by  simultaneously  rupturing  all  the  diaphragms.  The  total  volume  of  the  driver  system  will 
be  a  maximum  of  584  m3. 

Nitrogen  will  be  used  as  the  driver  gas.  The  nitrogen  is  to  be  stored  in  liquid  form  prior 
to  a  test.  When  needed,  the  liquid  nitrogen  will  be  pumped  into  a  pebble  bed  heater  2  where 
it  is  evaporated  and  heated  to  the  desired  temperature.  The  heated  nitrogen  is  then  sent  to 
the  driver  tubes.  The  current  gas  supply  design  allows  for  rapid  filling  of  the  driver  tubes 
to  minimize  loss  of  heat  to  the  uninsulated  driver  walls.  A  more  conventional  system  based 
on  compressors  would  fill  the  tubes  much  more  slowly.  To  avoid  heat  loss  from  the  gets  to 
the  driver  walls  with  a  compressor  based  system,  the  driver  walls  would  have  to  be  heated 
or  insulated,  neither  of  which  proved  to  be  practical. 

The  expansion  section  of  the  facility  will  be  semi-circular  with  a  cross-sectional  area  of 
163  m2  and  will  be  composed  of  reinforced  concrete  except  in  the  area  of  the  test  section, 
where  it  will  be  made  of  steel.  With  an  available  flow  area  of  this  size,  full-scale  military 
equipment  can  be  tested  under  realistic  nuclear  blast  and  thermal  conditions.  If  the  expan¬ 
sion  section  had  a  smaller  cross-sectional  area,  a  such  equipment  could  significantly  block 
the  flow  which  would  be  degrade  the  fidelity  of  the  blast  simulation. 

The  test  section  begins  89  m  from  the  beginning  of  the  expansion  section  and  is  17.7  m 
long.  This  location  of  the  test  section  is  considered  to  be  the  point  at  which  the  flow  from 
the  nine  driver  tubes  is  uniform  in  the  expansion  section.  A  thermal  radiation  source  will  be 
located  in  the  test  section  upstream  from  the  target.  The  thermal  radiation  associated  with 
a  nuclear  burst  is  simulated  through  the  combustion  of  a  mixture  of  powdered  aluminum 


and  oxygen.  The  combustion  products  are  removed  from  the  simulator  by  ejectors  or  vents 
located  above  the  thermal  radiation  source  in  the  test  section. 

An  active  Rarefaction  Wave  Eliminator  (RWE)  ^  will  be  located  at  the  downstream 
end  of  the  expansion  section.  This  device  is  used  to  eliminate  (or  minimize  the  effects 
of)  expansion  waves  which  are  generated  when  the  blast  wave  meets  the  open  end  of  the 
expansion  section.  These  unwanted  expansion  waves  travel  upstream  against  the  subsonic 
flow  and  toward  the  test  section.  They  destroy  the  fidelity  of  the  simulation  if  they  reach  the 
test  section.  The  RWE  is  termed  “active”  because  it  is  continuously  changing  the  available 
flow  area  so  that  the  expansion  waves  will  be  eliminated  for  the  entire  decaying  blast  wave. 
A  schematic  diagram  of  the  LB/TS  with  all  its  major  components  is  shown  in  Figure  1. 


II.  The  1:57  Scale  Large  Blast  Simulator 

Part  of  BRL’s  research  program  has  been  to  conduct  experiments  in  small  scale  shock 
tubes  in  an  effort  to  characterize  the  expected  gas  flow  in  the  LB/TS.  By  performing  exper¬ 
iments  in  a  small  scale  simulator,  a  large  number  of  inexpensive  tests  can  be  obtained  in  a 
relatively  short  period  of  time.  Also,  the  use  of  a  small  simulator  allows  for  easy  modification 
and  simple  operation. 

A  1:57  scale  LBS  (large  blast  simulator  with  no  thermal  radiation  capability)  has  been 
used  quite  extensively  at  BRL  to  study  the  expected  flow  in  the  LB/TS  and  the  effects  that 
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design  changes  will  have  on  the  flow.  Figure  2  is  an  illustration  of  the  small  scale  blast 
simulator.  Unlike  the  full  scale  LB/TS,  the  1:57  scale  LBS  has  a  single  driver  and  a  circular 
expansion  section.  It  models  the  nine  driver  LB/TS  by  lumping  the  available  flow  area  at 
any  given  location  into  a  single  diameter. 

The  model  was  configured  this  way  for  the  following  reasons: 

•  An  accurate  three  dimensional  model  would  be  difficult  to  construct  and  operate  because 
multiple  diaphragms  and  timing  for  diaphragm  breaking  would  be  difficult. 

•  A  three  dimensional  experimental  model  would  require  three  dimensional  computer 
simulations  which  are  complex  and  time  consuming. 

•  Pipe  sections  for  an  axisymmetric  model  are  readily  available. 

•  An  axisymmetric  model  can  be  simulated  with  one  and  two  dimensional  computer  codes 
which  require  fewer  computing  resources. 

With  this  two-dimensional  axisymmetric  model,  the  flow  characteristics  associated  with 
the  use  of  multiple  drivers  can  not  be  studied.  This  disadvantage,  however,  was  considered 
to  be  outweighed  by  the  advantages  of  the  simplified  system  described  earlier. 

A  rarefaction  wave  eliminator  was  not  used  during  the  tests  discussed  in  this  report. 
Rather,  a  very  long  expansion  section  was  employed.  The  use  of  a  long  expansion  section 
does  not  eliminate  the  rarefaction  wave.  It  simply  delays  the  arrival  of  the  rarefaction  wave 
at  the  test  section  until  after  the  period  of  interest  has  passed. 

Driver  gas  heating  was  employed  in  the  tests  but  the  method  by  which  the  gas  was 
heated  differed  from  that  intended  for  the  LB/TS.  Heating  of  the  driver  gas  was  accom¬ 
plished  by  placing  strip  heaters  on  the  outer  surface  of  the  driver  tube.  The  desired  driver 
gas  temperature  was  obtained  by  turning  the  strip  heaters  on  or  off  as  required.  One  purpose 
of  heating  the  driver  gas  is  to  provide  density  matching  on  either  side  of  the  contact  surface 
between  the  gas  originally  in  the  driver  and  the  shocked  gas  4.  The  other  purpose  of  driver 
gas  heating  is  to  add  energy  to  the  driver  gas  which  decreases  the  driver  pressure  required 
to  produce  a  given  shock  overpressure.  If  the  driver  gas  were  not  heated,  a  sudden  increase 
would  be  seen  in  the  dynamic  pressure  histories  as  the  contact  surface  passed  the  measure¬ 
ment  station.  Since  this  phenomenon  is  not  present  in  the  dynamic  pressure  histories  of  a 
free-field,  exponentially  decaying  blast  wave,  it  is  undesirable  in  the  LB/TS.  For  this  reason, 
an  effort  is  made  to  provide  contact  surface  density  matching  for  all  LBS  experiments. 


III.  Experimental  Measurements 


An  experimental  parametric  study  was  performed  in  the  1:57  scale  LBS  configured  as 
shown  in  Figure  2.  The  purpose  of  this  experimental  study  was  to  determine  the  perfor¬ 
mance  of  the  1:57  scale  LBS  in  simulating  free-field  nuclear  blast  waves  In  this  study, 
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the  driver  volume  and  initial  driver  pressure  were  varied  with  each  test.  In  all  the  tests, 
attempts  were  made  to  set  the  driver  temperature  according  to  the  contact  surface  density 
matching  criterion  discussed  in  Reference  4.  Table  1  lists  the  initial  test  conditions  for  all 
the  experiments  included  in  the  parametric  study. 


Table  1.  Test  Matrix 


Shot 

Driver 

Ambient 

Ambient 

Driver 

Driver 

Number 

Length 

Pressure 

Temperature 

Pressure 

Temperature 

(cm) 

(. kPa ) 

(A') 

(. kPa ) 

(AT) 

109 

87.31 

102.6 

296.05 

4378 

399.15 

110 

87.31 

102.3 

298.05 

6998 

453.15 

111 

87.31 

101.8 

297.05 

10170 

493.15 

112 

87.31 

102.0 

297.15 

12238 

499.15 

115 

87.31 

100.7 

295.85 

16030 

535.15 

116 

66.83 

102.3 

294.85 

4309 

393.15 

118 

66.83 

101.9 

295.65 

10617 

498.15 

119 

66.83 

101.8 

295.95 

13204 

520.15 

120 

66.83 

101.1 

296.15 

16651 

568.15 

126 

33.97 

102.2 

297.55 

4309 

387.15 

129 

33.97 

101.2 

296.05 

4412 

385.15 

127 

33.97 

101.2 

298.05 

6722 

449.15 

128 

33.97 

102.2 

298.55 

9308 

487.15 

130 

33.97 

101.4 

297.85 

11376 

510.15 

131 

18.97 

101.7 

297.75 

4344 

388.15 

132 

18.97 

101.5 

297.75 

6757 

460.15 

133 

18.97 

102.5 

297.05 

9997 

495.15 

134 

18.97 

103.1 

298.75 

11376 

503.15 

135 

18.97 

103.4 

297.15 

15686 

502.15 

136 

6.03 

102.6 

297.05 

4344 

386.15  | 

138 

6.03 

103.8 

297.65 

6895 

434.15 

137 

6.03 

102.8 

297.75 

9825 

489.15 

140 

6.03 

101.4 

297.35 

12928 

497.15 

139 

6.03 

103.3 

294.05 

15789 

530.15 

For  all  the  experiments  in  the  parametric  study,  pressure  measurements  were  made  at 
three  different  axial  locations  in  the  expansion  section.  Relative  to  the  upstream  end  of  the 
expansion  section,  the  locations  of  these  gages  are  76.2  cm,  127.0  cm  and  177.8  cm. 

Two  pressure  gages  were  located  at  each  of  the  three  axial  positions.  A  static  pressure 
gage  was  located  with  its  sensing  element  flush  with  the  interior  surface  of  the  expansion 
section  wall,  12.7  cm  from  the  axis  of  symmetry.  A  stagnation  probe  was  positioned  at  the 
half  radius  point,  6.35  cm  from  the  axis  of  symmetry,  and  was  facing  upstream  into  the  axial 
flow. 
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The  location  of  the  last  pair  of  gages  is  important  in  that  it  is  this  location,  when 
scaled,  which  corresponds  to  the  location  of  the  test  station  in  the  full  scale  LB/TS.  The 
data  recorded  at  this  station  was  ultimately  used  to  determine  the  shock  overpressure  and 
simulated  yield  for  each  test  in  the  experimental  parametric  study.  The  data  recorded  at 
this  station  is  the  most  significant  in  terms  of  LB/TS  modelling.  The  data  recorded  at  the 
remaining  stations  is  also  important  in  that  it  allows  the  analyst  to  see  how  the  flow  is 
changing  as  the  shock  wave  travels  down  the  expansion  section. 

By  measuring  the  stagnation  and  static  pressures  at  each  axial  location,  the  Mach 
number  and  dynamic  pressure  histories  can  be  obtained.  Accurate  dynamic  pressure  histories 
are  crucial  to  blast  simulation  because  it  is  the  dynamic  pressure  that  causes  targets  to 
experience  large  displacements  and  overturning.  In  fact,  the  intended  operating  envelope  of 
the  LB/TS  is  based  on  dynamic  pressure  impulse  for  shock  overpressures  above  70  kPa.  If 
the  shock  overpressure  is  below  this  value,  the  dynamic  pressure  history  becomes  difficult  to 
measure  experimentally. 

In  order  to  find  the  dynamic  pressure  one  must  first  determine  the  Mach  number,  which 
can  be  calculated  using  the  static  and  stagnation  pressures.  To  obtain  the  Mach  number 
history,  one  simply  repeats  the  following  set  of  calculations  for  each  point  in  the  pressure 
histories.  Equation  1  assumes  the  flow  is  subsonic  and  that  it  is  brought  to  rest  at  the 
probe’s  stagnation  point  through  an  isentropic  process  The  isentropic  relation  expressed 
in  Equation  1  can  be  solved  for  the  square  of  the  Mach  number  yielding  Equation  2. 

a-(i  + 

P  v 

■  (A) 

Where: 

•  M  is  the  local  Mach  number  of  the  gas, 

•  po  is  the  measured  stagnation  pressure, 

•  p  is  the  measured  static  pressure,  and 

•  7  is  the  ratio  of  specific  heats  (1.400). 

If  the  local  Mach  number  turns  out  to  be  greater  than  1.0  from  Equation  2,  then  the  flow 
is  considered  to  be  supersonic,  in  which  case  Equation  3  should  be  used.  Since  Equation  3 
has  no  closed  form  solution  for  A/,  a  Newton  iterative  solver  is  used  to  determine  its  value 
for  a  given  static  and  stagnation  pressure. 


( — r)/("r— l) 


(3) 
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Once  the  local  Mach  number  is  known,  the  dynamic  pressure,  q,  can  be  calculated  using 
Equation  4. 

q  =  'tpM2  (4) 

It  was  stated  earlier  that  for  low  shock  overpressures,  the  dynamic  pressure  becomes 
difficult  to  measure  experimentally.  As  the  shock  overpressure  is  reduced,  the  difference 
between  the  stagnation  pressure  and  the  static  pressure  becomes  small.  Thus,  the  ratio  of 
stagnation  pressure  to  static  pressure,  po/p,  approaches  unity.  As  illustrated  by  Equation  2, 
the  Mach  number  approaches  zero  as  the  pressure  ratio,  po/p,  approaches  unity.  Since 
the  dynamic  pressure  is  a  product  of  the  Mach  number  and  static  pressure,  as  defined  in 
Equation  4,  its  value  will  approach  zero  as  well.  This  phenomenon  is  illustrated  in  Figure  3, 
which  is  a  plot  of  stagnation  pressure  ratio,  po/poo,  and  dynamic  pressure  ratio,  q/poo,  as 
a  function  of  shock  pressure  ratio,  p/poo,  where  p^  is  equal  to  the  ambient  atmospheric 
pressure.  As  a  result  of  this  relationship  between  dynamic  pressure  and  shock  overpressure, 
small  disturbances  in  the  static  and  stagnation  pressure  histc'.ies  caused  by  electronic  noise 
or  turbulence  in  the  flow  will  be  exaggerated  in  the  experimental  dynamic  pressure  history. 

Effect  of  Shock  Overpressure  on  Dynamic  Pressure 


1.00  1.25  1.50  1.75  2.00  2.25  2.50  2.75  3.00  3.25  3.50 

Shock  Pressure  Ratio 


Figure  3.  Effect  of  Shock  Overpressure  on  Dynamic  Pressure  Ratio 
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IV.  Calculations 


To  complement  the  experimental  work,  a  computational  parametric  study  was  per¬ 
formed  using  the  BRL-Q1D  code  7.  The  BRL-QlD  code  was  chosen  for  its  relative  ease  of 
operation  and  quick  execution  time.  A  BRL-QlD  calculatioi  vas  performed  for  each  ex¬ 
periment  in  the  program.  This  computational  parametric  study  served  two  purposes.  First, 
by  comparing  the  experimental  results  to  the  calculated  results,  possible  anomalies  in  the 
experimental  records  could  be  detected  by  examining  those  cases  in  which  the  experimental 
and  calculational  data  suddenly  diverged.  Also,  the  experimental  data  could  be  used  to 
validate  the  code  so  that  further  calculations  could  be  performed.  This  validation  would 
make  it  possible  to  study  the  flow  characteristics  of  the  1:57  scale  LBS  for  cases  which  are 
beyond  the  operating  limits  of  the  actual  system. 

For  each  of  the  experimental  and  computational  results,  the  simulated  weapon  yield  was 
calculated.  This  was  done  by  multiplying  the  dynamic  pressure  impulse  by  the  scaling  factor 
of  57  and  then  applying  the  scaling  laws  presented  by  Glasstone  and  Dolan  ®.  These  scaling 
laws  relate  the  blast  parameters  between  different  explosive  yields.  The  resulting  explosive 
yield  is  expressed  in  kilotons  ( kT ).  By  comparing  the  experimental  and  computational  results 
to  a  known  standard  (in  this  case  1  KT)  of  equal  incident  overpressure,  the  simulated  yield 
for  the  test  can  be  determined  as  illustrated  in  Equation  5. 


h  \wj 


(5) 


In  this  equation,  I\  is  the  dynamic  pressure  impulse  for  the  reference  yield,  and  /2  is 
the  dynamic  pressure  impulse  at  the  same  pressure  level  for  the  predicted  yield,  W2-  Solving 
Equation  5  for  W2  yields  Equation  6.  Substitution  of  the  reference  yield  Wx  equal  to  1.0  AT 
leaves  the  predicted  yield  expressed  in  kilotons  as  a  function  of  the  measured  impulse  and 
the  impulse  of  the  reference  yield. 


(6) 


One  drawback  of  the  BRL-QlD  code  is  that  it  is  unable  to  accurately  predict  the  peak 
shock  overpressure  of  a  given  test.  It  has,  however,  demonstrated  the  ability  to  predict 
static  overpressure  impulse  and  dynamic  pressure  impulse  with  reasonable  accuracy.  The 
yield  calculations  discussed  above  are  derived  from  the  peak  shock  overpressure  and  the 
dynamic  pressure  impulse  of  a  particular  test.  Since  the  overpressure  results  obtained  from 
the  BRL-QlD  code  are  known  to  be  inaccurate,  the  simulated  weapon  yield  for  a  BRL-QlD 
calculation  is  computed  from  the  experimentally  measured  peak  shock  overpressure  and  the 
BRL-QlD  dynamic  pressure  impulse. 
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V.  Validation  of  the  BRL-Q1D  Code 


For  each  experiment  in  the  parametric  study,  a  matching  BRL-Q1D  calculation  was  per¬ 
formed.  A  simulated  explosive  yield  was  calculated  from  the  experimentally  measured  static 
overpressure  and  dynamic  pressure  impulse  for  the  experiments  and  calculations.  These  re¬ 
sults  are  summarized  in  Table  2.  The  second  column  in  this  table  lists  the  full  scale  driver 
volume  in  cubic  meters.  This  is  the  driver  volume  which  would  exist  if  the  simulator  shown 
in  Figure  2  were  enlarged  by  the  scaling  factor  of  57.  As  in  Table  1,  the  tests  are  listed  in 
groups  of  equal  driver  volume  with  the  largest  equal  to  1337  m3  and  the  smallest  of  117  m3. 
The  tests  in  each  group  are  listed  in  order  of  increasing  driver  pressure. 


Table  2.  Experimental  and  Calculational  Results 


Shot 

Number 

Full 

Scale 

Driver 

Volume 

(m3) 

Experimental 

Static 

Overpressure 

(kPa) 

Scaled 

Experimental 
Dynamic 
Pressure 
Impulse 
(kPa  -  s) 

Yield 
Based  on 
Experiment 
Results 

m 

Scaled 

Q1D 

Dynamic 
Pressure 
Impulse 
( kPa  —  s) 

Yield 
Based  on 
Q1D 
Results 

(*T) 

1337 

109 

1423.07 

HEBB 

'.;JS 

1337 

150 

2384.26 

IBfl 

111 

1337 

180 

■ 

7145.70 

n 

3195.42 

112 

1337 

205 

' 

7669.99 

N/A 

N/A 

115 

1337 

248 

91.200 

5320.72 

N/A 

N/A 

SB 

1030 

108 

12.825 

350.32 

9.063 

123.63 

UBS 

1030 

173 

33.060 

959.44 

42.522 

2051.88 

119 

1030 

215 

54.720 

1933.11 

63.156 

2972.09 

120 

1030 

250 

98.040 

6474.45 

N/A 

N/A 

126 

536 

N/A 

N/A 

4.560 

18.54 

129 

536 

4.805 

18.83 

127 

536 

138 

IH 

9.576 

55.46 

128 

536 

156 

14.364 

119.11 

16.587 

183.41 

130 

536 

178 

22.800 

282.38 

23.598 

313.08 

mm 

311 

103 

3.340 

7.51 

2.622 

3.63 

B 

311 

138 

6.498 

17.40 

5.301 

9.45 

Hi 

311 

166 

9.804 

29.77 

10.260 

34.12 

134 

311 

183 

11.400 

32.88 

12.215 

40.45 

135 

311 

223 

18.639 

69.18 

22.059 

114.67 

KB- 

117 

89 

1.254 

0.73 

0.798 

0.19 

IBM 

117 

110 

2.451 

2.37 

1.653 

0.73 

137 

117 

134 

2.850 

1.68 

2.907 

1.79 

140 

117 

152 

4.845 

4.95 

4.845 

4.95 

139 

117 

198 

5.586 

2.83 

6.544 

4.55 

The  data  in  Table  2  is  incomplete  for  tests  112,  115,  120  and  126.  In  tests  112,  115  and 
120,  the  computational  results  could  not  be  obtained  because  a  recompression  shock  located 
upstream  of  the  measurement  station  degraded  the  pressure  histories  at  the  measurement 
station.  This  made  the  recording  of  the  dynamic  pressure  impulse  for  these  calculations 
impossible.  The  effect  of  the  recompression  shock,  however,  was  not  observed  in  the  ex¬ 
perimental  records  for  these  tests.  During  test  126,  the  stagnation  pressure  gage  failed  to 
record  the  experimental  data.  This  made  the  calculation  of  the  dynamic  pressure  (and  hence 
dynamic  pressure  impulse)  impossible.  The  static  pressure  gage  functioned  properly  in  this 
test,  allowing  us  to  record  the  experimental  shock  overpressure  which  is  listed  in  Table  2  as 
104  kPa. 

For  each  of  the  five  driver  volumes  represented  in  Table  2,  the  results  are  illustrated 
in  a  chart  of  the  scaled  dynamic  pressure  impulse  versus  the  static  overpressure  for  the 
experimental  and  computational  results.  Figure  4  shows  the  results  of  the  tests  using  a  full 
scale  driver  volume  equal  to  1337  m3.  In  this  figure  the  dynamic  pressure  impulse  curve 
for  the  calculational  data  follows  that  of  the  experimental  data  but  at  a  lower  value.  As  a 
result  of  these  differences  in  impulse,  the  differences  in  the  yield  values  are  exaggerated  as 
evidenced  in  Table  2.  This  exaggeration  is  caused  by  the  nature  of  the  scaling  laws  described 
earlier. 

1:57  Scale  Heated  Driver  Experiments 
Full  Scale  Driver  Volume  =  1337  m’ 


Figure  4.  Dynamic  Pressure  Impulse  for  1337  m3  Driver  Volume 


As  illustrated  in  Equation  6,  the  simulated  yield  is  calculated  from  the  cube  of  the  ratio 
of  the  measured  impulse  and  the  reference  impulse.  Because  of  this  cubing  of  the  ratio, 
small  variations  in  the  measured  impulse,  I2,  may  lead  to  large  variations  in  the  calculated 
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yield,  In  addition,  as  shock  overpressure  decreases,  the  dynamic  pressure  impulse  for  the 
reference  yield  decreases.  As  a  result,  the  calculated  yield  is  more  sensitive  to  variations  in 
impulse  for  small  overpressures  than  for  large  overpressures.  Thus,  if  the  difference  between 
the  calculational  and  experimental  impulse  values  is  roughly  constant,  then  the  difference 
between  the  calculational  and  experimental  yields  will  be  largest  at  the  lowest  overpressures. 
For  these  reasons,  when  comparing  computational  data  to  experimental  data,  it  is  wise  to 
directly  compare  the  measured  impulse  rather  than  the  calculated  yield. 

The  results  of  the  tests  performed  with  an  equivalent  full  scale  driver  volume  of  1030  m3 
are  shown  in  Figure  5.  Once  again,  this  figure  compares  the  impulse  results  for  the  experi¬ 
ments  and  the  calculations.  As  in  the  previous  figure,  the  impulse  results  for  the  calculation 
seem  to  follow  the  trend  set  by  the  experiments.  The  differences  between  the  experimental 
and  computational  impulse  results  fall  between  15  and  29  percent. 

1:57  Scale  Heated  Driver  Experiments 
Full  Scale  Driver  Volume  =  1030  m1 


Figure  5.  Dynamic  Pressure  Impulse  for  1030  m3  Driver  Volume 


The  current  design  of  the  LB/TS  employs  a  maximum  driver  volume  of  584  m3.  The 
third  set  of  tests  in  the  parametric  study  was  performed  using  an  equivalent  full  scale  driver 
volume  of  536  m3.  The  results  from  this  set  of  tests  should  provide  a  good  estimate  of 
the  expected  performance  of  the  full  scale  LB/TS.  The  results  of  these  tests  are  shown  in 
Figure  6.  In  this  figure  one  can  see  that  the  dynamic  pressure  impulse  measured  in  the 
experiments  compares  very  well  with  the  computational  results.  At  138  kPa  this  error  is 
only  0.6  percent.  The  maximum  impulse  error  in  this  set  occurred  at  156  kPa  where  it  was 
15.5  percent. 

The  final  two  sets  of  tests  were  performed  using  equivalent  full  scale  driver  volumes  of 
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1 :57  Scale  Heated  Driver  Experiments 
Full  Scale  Driver  Volume  =  536  m’ 


Static  Overpressure  (kPa) 

Figure  6.  Dynamic  Pressure  Impulse  for  536  m3  Driver  Volume 

311  m3  and  117  m3.  The  results  of  these  tests  are  summarized  in  Figures  7  and  8.  For  the 
311  m3  tests,  the  impulse  error  was  a  minimum  of  4.6  percent  at  103  kPa  and  a  maximum  of 
21.5  percent  at  166  kPa.  The  results  for  the  117  m3  tests  can  be  seen  in  Figure  8.  In  these 
tests,  the  experimental  and  computational  data  compare  well  at  the  higher  overpressures 
where  the  minimum  error  is  0.0  percent  and  the  maximum  is  17.1  percent.  At  the  two  lowest 
overpressures,  however,  the  impulse  errors  are  greater  than  30  percent. 


VI.  Driver  Temperature  Requirements 

It  was  mentioned  earlier  that,  in  performing  the  experiments,  an  effort  was  made  to 
conduct  the  tests  in  accordance  with  the  contact  surface  density  matching  relationship  de¬ 
scribed  by  Opalka  in  Reference  4.  The  ability  of  the  experimentalist  to  meet  this  criterion 
will  now  be  examined. 

The  density  matching  criterion  is  described  in  terms  of  driver  pressure  ratio  and  driver 
temperature  ratio.  The  driver  pressure  ratio  is  defined  as  the  sum  of  the  driver  overpressure 
plus  the  ambient  pressure  divided  by  the  ambient  pressure.  The  driver  temperature  ratio  is 
defined  as  the  absolute  driver  gas  temperature  divided  by  the  absolute  ambient  temperature. 
The  relationship  between  the  driver  pressure  ratio  and  the  driver  temperature  ratio  for  the 
test  conditions  is  available  in  Table  3.  In  this  table  one  can  see  that  for  many  of  the 
tests,  the  driver  temperature  ratio  measured  from  the  experiments  was  significantly  less 
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Scaled  Dynamic  Pressure  Impulse  (kPa-s) 


1 :57  Scale  Heated  Driver  Experiments 
Full  Scale  Driver  Volume  =  311  m’ 


t - 1 - 1 - 1  i  i 

100  125  150  175  200  225 

Static  Overpressure  (kPa) 


Figure  7.  Dynamic  Pressure  Impulse  for  311  m3  Driver  Volume 


1:57  Scale  Heated  Driver  Experiments 
Full  Scale  Driver  Volume  =  117  m* 


75  100  125  150  175  200  225 

Static  Overpressure  (kPa) 


Figure  8.  Dynamic  Pressure  Impulse  for  117  m 3  Driver  Volume 
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than  the  ideal  temperature  ratio  as  defined  in  Reference  4.  This  difference  between  the 
experimental  and  ideal  temperature  ratios  is  illustrated  in  Figure  9,  which  is  a  chart  of 
the  driver  temperature  ratio  versus  the  driver  pressure  ratio  for  the  experimental  and  ideal 
conditions.  This  figure  illustrates  that  the  difference  between  the  experimental  and  ideal 
temperature  ratio  is  greatest  for  the  largest  driver  pressure  ratios.  From  this  chart,  one  can 
conclude  that  the  heating  apparatus  used  in  the  experiments  was  insufficient  for  heating  the 
driver  gas  to  the  highest  desired  temperatures. 


Table  3.  Experimental  and  Ideal  Driver  Temperature  Ratios 


Shot 

Number 


Driver 

Pressure 

Ratio 


43.6716 

69.4066 

100.9018 

120.9804 

160.1857 


43.1212 

105.1904 

130.7053 

165.69S3 


43.1624 

44.5968 

67.4229 

92.0763 

113.1893 


43.7139 

67.5714 

98.5317 

111.3395 

152.7021 


43.3392 

67.4258 

96.5739 

128.4951 

153.8461 


Experimental 

Driver 

Temperature 

Ratio 


1.3483 

1.5204 

1.6602 

1.6798 

1.8089 


.3334 
.6849 
1.7576 
1.9185 


.3011 

.3010 

.5070 

.6317 

.7128 


1.3036 

1.5454 

1.6669 

1.6842 

1.6899 


.2999 
.4586 
1.6428 
1.6719 
1.7999 


Ideal 

Driver 

Temperature 

Ratio 


1.3036 

1.6074 

2.0107 

2.1888 

2.4475 


1.2965 

2.0565 

2.2530 

2.4839 


1.2965 

1.3109 

1.5810 

1.9038 

2.1319 


1.3036 

1.5810 

1.9822 

2.1167 

2.3981 


1.3013 

1.5810 

1.9587 

2.2384 

2.4057 


The  primary  purpose  of  heating  the  driver  gas  is  to  provide  density  matching  on  either 
side  of  the  contact  surface  between  the  gas  originally  in  the  driver  and  the  shocked  gas. 
The  secondary  purpose  is  to  decrease  the  driver  pressure  required  to  produce  a  given  shock 
overpressure.  If  the  driver  gas  is  not  sufficiently  heated,  a  sudden  increase  will  be  present 
in  the  dynamic  pressure  histories.  This  effect  can  be  seen  in  Figure  10  which  is  a  dynamic 
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1 :57  Scale  Heated  Driver  Experiments 
Experimental  and  Ideal  Driver  Temperature  Ratios 


Figure  9.  Experimental  and  Ideal  Driver  Temperature  Ratios 

pressure  history  from  a  BRL-Q1D  calculation  using  the  experimental  driver  temperature 
ratio.  In  this  figure,  the  peak  dynamic  pressure  is  about  120  kPa  and  the  dynamic  pressure 
impulse  is  about  0.70  kPa  —  s.  The  sudden  increase  in  dynamic  pressure  is  labeled  as  a  “Cold 
Gas  Effect”  in  the  figure.  This  effect  is  caused  by  the  passing  of  the  contact  surface  through 
the  measurement  station.  This  cold  gas  effect  creates  a  net  increase  in  dynamic  pressure 
impulse. 

Another  dynamic  pressure  history  is  shown  in  Figure  11.  This  history  is  the  result  of  a 
subsequent  BRL-QlD  calculation  which  employed  the  ideal  driver  temperature  ratio.  As  a 
direct  result  of  the  additional  driver  gas  heating,  the  peak  dynamic  pressure  has  increased  to 
about  140  kPa  and  the  cold  gas  effect  is  gone.  The  dynamic  pressure  impulse  for  this  calcu¬ 
lation  is  about  0.66  kPa  —  s  which  is  less  than  that  of  the  calculation  using  the  experimental 
temperature  ratio. 

Since  the  additional  heating  caused  the  incident  shock  overpressure  to  increase  and  the 
dynamic  pressure  impulse  to  decrease,  the  simulated  yield  for  the  calculation  using  the  ideal 
temperature  ratio  will  be  less  than  that  which  used  the  experimental  temperature  ratio. 
For  this  reason,  performance  predictions  for  the  full  scale  LB/TS  which  are  derived  from  the 
experimental  results  discussed  earlier  may  not  be  accurate  for  the  higher  shock  overpressures. 
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Time  (ms) 

Figure  10.  Dynamic  Pressure  Using  Experimental  Temperature  Ratio 


VII.  Calculations  Using  Ideal  Driver  Temperatures 


In  an  effort  to  develop  a  better  estimate  of  the  expected  performance  of  the  LB/TS, 
a  second  computational  parametric  study  was  performed  with  the  BRL-QlD  code.  This 
parametric  study  employed  the  ideal  driver  temperature  ratio  for  each  test  rather  than  the 
experimental  driver  temperature  ratio.  No  additional  calculations  were  performed  for  the 
low  overpressure  cases  in  which  the  experimental  and  ideal  driver  temperature  ratios  were 
very  close. 

The  results  of  this  second  computational  parametric  study  are  summarized  in  Table  4. 
In  this  table,  the  tests  for  which  a  second  calculation  was  performed  are  noted  with  an 
asterisk  (*)  next  to  the  shot  number.  As  in  Table  2,  the  computational  data  for  test  115  was 
not  obtained  because  of  the  presence  of  a  recompression  shock  which  destroyed  the  pressure 
histories  at  the  measurement  station.  However,  the  data  for  tests  112  and  120  were  obtained 
in  the  second  computation  study,  whereas  they  were  lost  in  the  first  study  due  to  the  effect 
of  a  recompression  shock.  Apparently  the  additional  driver  gas  heating  used  in  the  second 
calculation  helped  to  suppress  the  effect  of  this  recompression  shock  at  the  measurement 
station.  This  phenomenon  is  documented  in  Reference  1. 

Two  figures  are  given  for  each  set  of  tests  corresponding  to  a  given  driver  volume.  The 
first  of  each  pair  of  figures  is  a  plot  of  the  scaled  dynamic  pressure  impulse  versus  the  static 
overpressure.  The  second  figure  in  each  pair  is  a  chart  of  the  simulated  weapon  yield  versus 
the  static  overpressure.  These  charts  are  presented  in  Figures  12  through  21. 

Using  the  ideal  driver  temperature  ratio  in  these  calculations  eliminated  the  cold  gas 
effect  discussed  earlier.  As  a  result,  the  dynamic  pressure  impulse  for  these  calculations  is 
lower  than  that  obtained  from  the  previous  set  of  calculations.  This  reduction  in  dynamic 
pressure  impulse  is  evident  in  the  charts  which  plot  the  impulse  against  shock  overpressure 
(Figures  12,  14,  16,  18  and  20).  It  is  interesting  to  note  that  the  relationship  between 
dynamic  pressure  impulse  and  shock  overpressure  is  almost  linear  for  all  five  sets  of  tests. 

It  was  previously  stated  that  heating  the  driver  gas  to  the  ideal  temperature  ratio  would 
cause  the  initial  yield  estimates  to  decrease.  As  illustrated  in  Figures  13,  15,  17,  19  and  21 
the  predictions  of  simulated  weapon  yield  from  the  second  study  are  significantly  less  than 
those  of  the  earlier  study.  For  example,  the  yield  for  test  119  drops  from  2972  kT  in  the 
earlier  study  to  1288  k T  in  the  second  study.  This  reduction  in  yield  is  evident  in  other  tests 
as  well.  It  is  also  interesting  to  note  that  in  all  the  yield  versus  shock  overpressure  figures, 
the  curves  are  approximately  asymptotic  in  nature  (their  slopes  decrease  with  increasing 
shock  overpressure). 
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Table  4.  Results  of  Q1D  Calculations  Using  Ideal  Driver  Temperature 


Shot 

Number 

Full 

Scale 

Driver 

Volume 

(m3) 

Expected 

Static 

Overpressure 
( kPa ) 

Scaled 
QlD 
Dynamic 
Pressure 
Impulse 
( kPa  —  s) 

Yield 
Based  on 
QlD 
Results 

m 

■n 

| 

■ 

j 

■111111 

in* 

1337 

1801.50 

112* 

' 

227 

60.990 

2243.21 

115* 

Km 

N/A 

N/A 

N/A 

116 

108 

9.063 

123.63 

118* 

184 

38.190 

1161.93 

119* 

1030 

235 

53.010 

1288.04 

1030 

275 

75.240 

2094.39 

536 

4.560 

18.54 

■SrSS 

536 

■ES 

4.805 

18.83 

127 

536 

138 

9.576 

55.46 

128* 

536 

163 

16.530 

154.48 

130* 

536 

192 

22.230 

195.41 

WtSM 

mm 

103 

2.620 

3.63 

K9 

m 

138 

5.301 

9.45 

133* 

311 
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*  Calculations  performed  in  second  QlD  study. 
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Scoled  Dynomlc  Pressure  Impulse  (kPa-s) 


Q1D  Results  Using  Ideal  Temperature  Ratio 
Full  Scale  Driver  Volume  =  1337  m* 


Figure  12.  Impulse  for  1337  m3  Driver  Using  Ideal  Temperature  Ratio 


Q1D  Results  Using  Ideal  Temperature  Ratio 
Full  Scale  Driver  Volume  =  1337  m* 
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Simulated  Weapon  Yield  (kT) 


Q1D  Results  Using  Ideal  Temperature  Ratio 


Static  Overpressure  (kPa) 
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mic  Pressure  Impulse  (kPa-s) 


Q1D  Results  Using  Ideal  Temperature  Ratio 
Full  Scale  Driver  Volume  =  536  m* 


Figure  16.  Impulse  for  536  m3  Driver  Using  Ideal  Temperature  Ratio 


Q1D  Results  Using  Ideal  Temperature  Ratio 
Full  Scale  Driver  Volume  =  536  m* 


Scaled  Dynamic  Pressure  Impulse  (kPa- 


01 D  Results  Using  Ideal  Temperature  Ratio 
Full  Scale  Driver  Volume  =  311  m* 


Figure  18.  Impulse  for  311  m3  Driver  Using  Ideal  Temperature  Ratio 


01 D  Results  Using  Ideal  Temperature  Ratio 
Full  Scale  Driver  Volume  =  311  m* 
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Scaled  Dynamic  Pressure  Impulse  (kPa-s) 


VIII.  Conclusion 


This  report  has  attempted  to  analyze  the  results  of  two  calculational  studies  performed 
using  the  BRL-Q1D  code.  The  first  study  attempted  to  match  the  results  obtained  from  am 
experimental  parametric  study.  The  purpose  of  this  experimental  study  was  to  validate  the 
ability  of  the  BRL-Q1D  code  to  predict  the  static  and  dynamic  pressure  impulse  in  the  1:57 
scale  Large  Blast  Simulator.  The  results  of  the  first  computational  parametric  study  were 
analyzed  and  found  to  be  in  reasonable  agreement  with  the  experimental  results. 

During  the  experimentation,  attempts  were  made  to  heat  the  driver  gas  according  to  a 
predetermined  relationship  between  driver  gas  pressure  and  driver  gas  temperature.  When 
comparing  the  experimental  conditions  to  the  ideal  pressure/temperature  relationship,  the 
driver  gas  was  found  to  have  been  insufficiently  heated  for  many  of  the  high  pressure  tests. 
Failure  to  fully  heat  the  driver  gas  can  artificially  increase  the  measured  dynamic  pressure 
impulse.  This  effect  is  subtle  and  can  go  unnoticed  in  experimentally  measured  dynamic 
pressure  records. 

The  nature  of  the  scaling  laws  which  are  used  to  calculate  the  simulated  weapon  yield 
indicates  that  tests  with  insufficient  driver  gas  heating  may  produce  overly  optimistic  esti¬ 
mates  for  the  performance  of  the  LB/TS.  In  an  attempt  to  improve  the  performance  esti¬ 
mates,  a  second  computational  parametric  study  was  performed.  This  study  employed  the 
ideal  driver  gas  pressure/temperature  relationship  in  the  initial  conditions.  The  results  of 
this  study  prove  that  the  initial  performance  estimates  were  indeed  optimistic  due  to  the 
insufficient  driver  gas  heating.  The  revised  performance  predictions  which  resulted  from 
this  second  study  were  significantly  lower  than  those  from  the  experiments  and  the  initial 
computational  parametric  study. 
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San  Leandro,  CA  94577 

2  R&D  Associates 

ATTN:  Technical  Library 
Dr.  Allan  Kuhl 
P.O.  Box  9695 
Marina  Del  Rey,  CA  90291 

1  R&D  Associates 
ATTN:  G.P.  Ganong 
P.O.  Box  9377 
Albuquerque,  NM  87119 

2  Science  Applications,  Inc. 

ATTN:  W.  Layson 

John  Cockayne 
P.O.  Box  1303 
1710  Goodridge  Drive 
McLean,  VA  22102 

1  Science  Applications  International 

Corp. 

ATTN:  Mr.  J.  Guest 
2301  Yale  Blvd.  SE  Suite  E 
Albuquerque,  NM  87106 

1  Sparta,  Inc. 

Los  Angeles  Operations 
ATTN:  I.  B.  Osofsky 
3440  Carson  Street 
Torrance,  CA  90503 

1  Sunburst  Recovery,  Inc. 

ATTN:  Dr.  C.  Young 
P.O.  Box  2129 

Steamboat  Springs,  CO  80477 

1  Sverdrup  Technology,  Inc. 

ATTN:  R.  F.  Starr 
P.O.  Box  884 
Tullahoma,  TN  37388 


3  SRI  International 

ATTN:  Dr.  G.  R.  Abrahamson 
Dr.  J.  Gran 
Dr.  B.  Holmes 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

2  S-CUBED 

A  Division  of  Maxwell 
Laboratories,  Inc. 

ATTN:  C.  E.  Needham 

Dr.  Lynn  Kennedy 
2501  Yale  Blvd.,  SE 
Albuquerque,  NM  87106 

3  S-CUBED 

A  Division  of  Maxwell 
Laboratories,  Inc. 

ATTN :  Technical  Library 
R.  Duff 
K.  Pyatt 
PO  Box  1620 
La  Jolla,  CA  92037-1620 

1  Texas  Engineering  Experiment 

Station 

ATTN:  Dr.  D.  Anderson 

301  Engineering  Research  Center 

College  Station,  TX  77843 

1  Thermal  Science,  Inc. 

ATTN:  R.  Feldman 
2200  Cassens  Dr. 

St.  Louis,  MO  63026 

1  TRW 

Ballistic  Missile  Division 
ATTN:  H.  Konnan, 

Mail  Station  526/614 
P.O.  Box  1310 
San  Bernadino,  CA  92402 
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1  Battelle 

ATTN:  TACTEC  Library, 

J.  N.  Higgins 
505  King  Avenue 
Columbus,  OH  43201 

1  California  Institute  of  Technology 
ATTN:  T.  J.  Ahrens 

1201  E.  California  Blvd. 

Pasadena,  CA  91109 

2  Denver  Research  Institute 
ATTN:  Mr.  J.  Wisotski 

Technical  Library 
P.O.  Box  10758 
Denver,  CO  80210 

1  Massachusetts  Institute  of  Technology 

Aeroelastic  and  Structures  Research 
Laboratory 

ATT*  :  Dr.  E.  A.  Witmer 
Cambridge,  MA  02139 

1  Massachusetts  Institute  of  Technology 

ATTN:  Technical  Library 
Cambridge,  MA  02139 


4  Southwest  Research  Institute 

ATTN:  Dr.  W.  E.  Baker 
A.  B.  Wenzel 
Dr.  C.  Anderson 
S.  Mullin 

6220  Culebra  Road 
San  Antonio,  TX  78284 

1  Stanford  University 

ATTN:  Dr.  D.  Bershader 
Durand  Laboratory 
Stanford,  CA  94305 

3  University  of  Maryland 

Department  of  Mechanical  Engineering 
ATTN:  Dr.  W.  Fourney 
Dr.  R.  Dick 
Dr.  J.  Williams 
College  Park,  MD  20742 

Aberdeen  Proving  Ground 
1  Cdr,  USATECOM 

ATTN:  AMSTE-TE-F,  L.  Teletski 

1  Cdr,  USATHMA 

ATTN:  AMXTH-TE 


2  New  Mexico  Engineering  Research 
Institute  (CERF) 

University  of  New  Mexico 
ATTN:  Dr.  J.  Leigh 
Dr.  R.  Newell 
P.O.  Box  25 

Albuquerque,  NM  87131 


1  Northrop  University 

ATTN:  Dr.  F.  B.  Safford 
5800  W.  Arbor  Vitae  St. 
Los  Angeles,  CA  90045 
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